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I. INTRODUCTION 



An ionization process is removal of one or more electrons from a quantum system, such 
as atom, molecule, quantum heterostructure etc., in any manners. If the ionization pro- 
cess is induced by a photon, this process is called photoionization. The photoionization of 
atomic systems provides the opportunity to investigate the dynamic interplay among many 
body electron-electron correlations and relativistic effects.- It is the simplest process giving 
detailed information on the atomic and molecular structure.- Numerous experimental and 
theoretical studies on photoionization cross section of different atoms have been published by 
many authors.- - - Photoionization cross sections are of great significance due to its numerous 
applications in space research, astrophysics, radiation protection, laser designing, controlled 
thermonuclear research and different types of laboratory plasma.-^ 

On the other hand, because quantum dot (QD) heterostructures exhibit atomic properties 
such as energy levels, density of states etc., they are called as artificial atoms. 6 The investi- 
gation of photoionization cross section of a hydrogenic impurity in a QD is very important 
especially for different optoelectronic device applications such as quantum dot infrared pho- 
todetectors. Although there are no experimental studies, the photoionization cross section 
of a hydrogenic impurity in low- dimensional quantum structures has been studied theo- 
retically by many authors.- - — The cross section strongly depends on the impurity binding 
energy and its wave function.- The photoionization cross section is still studied extensively 
under different physical parameters (i.e electric field, magnetic field, hydrostatic pressure) 
in low-dimensional quantum nano structures. Recently, §ahm£ has investigated the pho- 
toionization cross section of both D° and D~ impurities in a QD for an on-center impurity 
as a function of the dot radii and the normalized photon energies. He demonstrated that 
the cross section drastically depends on the dot size, the photon energy and the number 
of electrons in the QD. Barseghyan et al.— studied the behavior of the binding energy and 
the photoionization cross-section of a donor impurity in a cylindrical-shape GaAs/AlGaAs 
quantum dot, under the effects of hydrostatic pressure and applied electric and magnetic 
fields in-growth direction. They showed that when the hydrostatic pressure is increased, 
the photoionization cross section grows. Ham and Spector- analyzed the dependence of the 
photoionization cross section on the energy and polarization of the photons in a spherical 
QD as a function of dot radius and impurity locations for both infinite and finite potential 
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barriers. They found that the cross section is independent of the polarization of the photons 
for an on-center impurity, while it depends on the polarization of the photon field for an 
off-center one. Sajeev and Moiseyev 15 investigated the autoionization and photoionization 
in two-electron spherical quantum dots. Baldea and Koppel 16 showed that the photoioniza- 
tion is a valuable method, which can be used along with transport measurements, to study 
single electron transistors. Most recently, Lin and Ho^ have studied the photoionization 
cross section of a hydrogen impurity in spherical quantum dots which are modelled by fi- 
nite oscillator and Gaussian potentials. All of these studies have been reported for a single 
core/shell QD. 

The production of the spherical multi-layered quantum dot (MLQD) has become possible 
using wet chemical techniques.— 1 ^ The electronic properties of multilayered spherical quan- 
tum dots under various physical effects have been investigated by some authors.— ~— These 
studies are related to electronic energy level and impurity binding energy calculations. It is 
seen from these studies that the energy levels and impurity binding energies are very sensi- 
tive to the layer thickness. Until now, to the best of our knowledge, although the electronic 
structure has been studied by many authors, the photoionization cross section properties 
have not been investigated in a MLQD. The investigation of photoionization cross section in 
MLQD may give rise to an opportunity for producing a new generation of quantum devices 
such as laser source, photodetectors, single electron transistor etc. For example, a single 
electron transistor can work as depending on the photoionization cross section.— The con- 
trollability of the shell thicknesses may provide different advantages in the production of 
these kind of devices working based on the photoionization cross section. 

The main aim of this work is to calculate the photoionization cross section of a hydrogenic 
impurity in a spherical multi-layered semiconductor quantum dot, and to investigate the de- 
pendence of the photoionization on the layer thickness. For this purpose, the ground and 
first excited energy levels and corresponding wave functions are calculated by solving the 
Schrodinger equation with shooting method for an electron confined in a MLQD. The vari- 
ations of the impurity binding energies and the photoionization cross section are examined 
as a function of shell thicknesses and photon energies. 

The rest of this paper is organized as follows: In the next section, the model and calcu- 
lations are presented. The calculation results and discussion of their probable reasons are 
given in section III. A brief conclusion is presented in the last section. 
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II. MODEL AND CALCULATION 



In this study, a spherically symmetric core/shell/ well/shell quantum dot structure is 
considered. As seen from Fig. (TJ GaAs core is a sphere of radius R±. The well region is 
that another GaAs spherical shell of inner radius R 2 and outer radius -R3 is concentric with 
the core region. The gap between the core and well regions is filled with AlGaAs as a shell 
(barrier) zone which isolates the core and well regions and also allows electron to tunnel 
between these regions. The shell thickness is T s = R 2 — R\ and the well region thickness is 
T w = R% — R 2 . All structure is embedded in an AlGaAs matrix material. 

In the effective mass approximation, for a spherically symmetric quantum dot, the single- 
particle Schrodinger equation is given as 



rrr r 



Ze 2 



H — h V(r) 

2m* (r) r 2 



Rn,e (r) = e n , e R n/ (r) , (1) 



where h is reduced Planck constant, m*(r) is the position-dependent electron's effective 
mass, Z is charge of the impurity, n(r) is the position dependent dielectric constant, r*j is 
the position of the impurity, i is the angular momentum quantum number, V(r) is the finite 
confining potential, e n ^ is the electron energy eigenvalue, and R n ^(r) is the radial wave 
functions of the electron. It should be noted that, in the case without the impurity, the 
charge of the impurity is Z = and in the existence of the impurity, it is Z — 1. In the 
calculations, we have taken into consideration an on-center impurity and hence the r*j = 0. 
The mathematical expression of the confining potential is 



V(r) 



0. 



r < R x and R 2 < r < R 3 



(2) 



V b , Ri<r < R 2 and r > R 3 

In order to determine the single particle energy levels and corresponding wave functions, 
Eq.(Tjrj) is full numerically solved by using the shooting method. As is well known, this 
technique converts an eigenvalue problem to an initial value problem. For this purpose, 
Hamiltonian operator is discretized on a uniform radial mesh in ID using the finite differ- 
ences, then Eq.flT]) can be reduced to an initial value equation. Here, Ar width between two 



mesh points is chosen as 0.005. All details of this method can be found in Ref. 
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FIG. 1. Schematic representation and potential profile of a core/shell/well/shell QD. 



A. Theory of photoionization cross-section 

The photoionization process is an optical transition that takes place from the impurity 
ground state as the initial state to the conduction subbands (continuum), which starts 
above the confining potential of the dot. The sufficient energy is required in order for the 
transition to occur-^. The photoionization cross section can be defined as the ionization 
probability of the electrons from the bound state under an external optical excitation. It 
is strongly dependent on the confinement potential and the photon energy. The excitation 
energy depends on the photoionization cross section associated with an impurity and, as 
in the bulk case, it can be calculated in frame of the Fermi's golden rule and the dipole 
approximation as following 24 



where n r is the refractive index of the semiconductor, k is the dielectric constant of the 
medium, (3fs = e 2 / He is the fine structure constant, and Hcj is the photon energy. F e ff/F 
is the ratio of the effective electric field F e ff of the incoming photon and average field F 
in the medium 2 ^. < ipi\r\ipf > is the matrix element between the initial and final states of 
the dipole moment of the impurity, ipf and ipi are the wave functions of the final and initial 
states, and Ef and E { are the corresponding energy eigenvalues of these states, respectively. 
These eigenvalues and eigenvectors are determined by the shooting method. 

In a spherical QD, the selection rules (A£ = ±1) determine the final state of the electron 
after the impurity photoionization. The Is impurity state is taken as ground level for D° 
and the lp level of one electron system without the impurity is taken as the final state for 
D° as similar with previous studies.- 1 ^— ^ 



2 n r An 2 



f3 FS hu Mi kl if>f)\ 2 S (E f -Ei-hJ), 



(3) 
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In order to calculate the numerical values of the photoionization cross-section given by 
Eq.fllJ), F e ff/F is taken as approximately unity because the calculation of this quantity is 
very difficult and it has no effect on the photoionization cross-section shaped— . Also, the 
initial and final state wave functions (ipi and ipf) are determined by the multiplication of 
the radial wave function with the spherical harmonics i.e. R n ^(r)Y^^ m {9 , if). Here, Y^ m (#, ip) 
is the spherical harmonics. In addition, the 5-function in Eq.Q is replaced by a narrow 
Lorentzian by means of, 



S (E f -Ei-fvjj) 



hv 



vr [(Huj - (E f - Ei)) + (hT) 2 ] 
Here, T is the hydrogenic impurity linewidth and taken as 0.1R*. 



(4) 



III. RESULTS AND DISCUSSION 

The atomic units have been used throughout the calculations, where h = m = e = 1. 
Effective Bohr radius is ~ 100 A and effective Rydberg energy is R* ~ 5.25 meV. The 
confinement potential has been taken as V& = 228 meV. The material parameters have 
been taken as m GaAs = 0.067 m Q , m A iGaAs = 0.088 m , K GaAs = 13.18, K A i GaAs = 12.8. 
Also the effective masses of electrons inside GaAs and AlGaAs are m* and m^, and the 
dielectric constants are K\ and k 2 , respectively. The position-dependent effective mass and 
the dielectric constant may be defined as follows 26 



1, r <Ri and R 2 < r < R 3 



m*(r) = < 



Ri < r < R 2 and r > R 3 



K\T ) = < 



1, r <R l and R 2 < r < R 3 



f 2 , Ri <r < R 2 and r > R 3 



(5) 



The binding energy of a neutral donor impurity [D°) in a QD is given as 6 



E b (D°)=E 10 -E 10 (D°). 



(6) 
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FIG. 2. Impurity binding energy of a multi-layered QD as a function of core radius, R±, for constant 
T s and T w . 

Here, E w is ground state energy of the single electron QD without the impurity and £'io(-D ) 
is that of the neutral donor impurity. 

Figure [2] shows the variation of the binding energy of the D° impurity as a function of 
the core radius of the MLQD, Ri. Here, the shell thickness is T s = 0.4 and the well width 
is T w = 0.4 <2q. As seen from the figure, in small core radius, the binding energy increases 
rapidly with increasing R±. It reaches a maximum value at about Ri = 0.6 and decreases 
smoothly with further increasing of the R\ and goes to a constant value. This behavior 
is rather similar to the binding energy of a neutral hydrogenic donor impurity in a single 



Figure [3] shows the binding energy of the impurity as a function of the shell thickness 
for T w = 0.4 Oq, and two different core radii, R± = 0.3 and 0.6 a^. As seen from the 
figure, the behaviors of the binding energies are completely different from the each other. 
When the core radius is Ri = 0.3 a^, the binding energy decreases with increasing shell 
thickness. On the other hand, when the core radius is Ri = 0.6 Oq, the binding energy 
increases firstly and remains at a constant value for almost T s > 0.6 a^. These changes are 
explained as follows: In the case of small core radius, the electron is confined in the well 
region for both cases without [Z = 0) and with {Z = 1) the impurity because the core 
region is rather small. However, in the presence of the impurity, the electron feels the effect 
of the impurity at smaller shell thickness and hence the energy level is pulled down because 
of the attractive coulomb potential of the impurity. As a result of this, the binding energy 
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FIG. 3. (Color online) Impurity binding energy of a core/shell/well/shell QD as a function of shell 
thickness for constant well width and two different core radii. 

becomes larger at small shell thickness. In contrast, while the shell thickness increases, 
the attractive effect of the impurity on the electron decreases so that, the electron is less 
affected from the impurity. These behaviors are clearly seen from the left panel of FigJU 
Accordingly, the binding energy decreases with increasing T s . By further increase in the T s , 
the electron is confined completely in the well region for both Z = and Z = 1 cases, and 
it is insubstantially affected from the impurity. In the case of large core radius, while the 
electron is confined in both the core and well regions for Z = case, it is almost confined in 
the core region for Z = 1 case. As a result, the energy difference is relatively small. When 
the shell thickness, T s , is larger, the electron is completely confined in the core region for 
both Z = and Z = 1 cases and the effect of the well region vanishes as seen from the 
right panel of FigJU Hence, the attractive coulomb effect of the impurity becomes more 
significant which results in the larger binding energy as seen in Fig. [3j 

We expect that the impurity binding energy will be affected drastically from the width 
of well region. In Fig. [5j the impurity binding energies are plotted as a function of the well 
width, T w , for two different shell thickness at R\ = 0.8 a^. When we look at the figure, we 
see similar treatments in the binding energies for both shell thicknesses. Nevertheless, we 
observe that the variation of the binding energy is smoother for T s = 0.2 as compared 
with case of T s = 0.4 a^. The physical reasons of these treatments can be explained as 
follows: For the first value of the shell thickness, in spite of the fact that the electron is 
confined inside the core region for both Z = and Z = 1 cases, the well region influences 
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FIG. 4. (Color online) Ground state energies and probability densities at same T w = 0.4 for 
Z = and Z = 1 cases and different layer thickness. The other layer parameters are given on the 
figure. 
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FIG. 5. (Color online) Variation of the impurity binding energy as a function of the well widths, 
T w , for R x = 0.8 ag, and T s = 0.2 and 0.4 og. 



the energy levels even for small values of the well width. Moreover, the energy difference 
between Z = and Z = 1 cases is larger, so that the binding energy increases. When 
the well region becomes much larger, the electron tunnels easily to the well region and the 
energy levels approaches to each other. This situation is shown in left panel of FiglHJ On 
the other hand, when the shell thickness is T s = 0.4 <Zq, the tunneling is not so easy and the 
electron is confined almost completely in the core region in comparison with the T s = 0.2 Oq 

9 




FIG. 6. (Color online) Ground state energies and probability densities at same Ri = 0.8 Oq for 
Z = and Z = 1 cases and different layer thickness. 



case as seen in right panel of FigJBl As opposed to the T s = 0.2 czq case, when the well 
region is much larger, the electron tunnels completely to the well region and as depending 
on this, the decreasing in the binding energy becomes sharper with increasing shell thickness 
as observed in Fig. [5j 

The photoionization cross section of a donor impurity is an important phenomena espe- 
cially for device applications of quantum heterostructures. If the photon energy is equal 
to the energy difference between the initial and the final states (i.e. resonant energy, 
fauu = Ef — Ei), all terms of EqOJ except that photon energy and dipole matrix element, con- 
tribute as a constant to the photoionization cross section. In this resonant photon energy 
case, the Lorentz function becomes also another constant as 1/nhT. On the other hand, 
the energy difference between the initial and the final states and the overlapping of their 
wave functions (hence the dipole matrix element) play an effective role on the shape of the 
photoionization cross section.- 

Fig. [7] shows the variation of the peak values of the photoionization cross section with 
the core radius, R\, for constant shell and well widths, i.e. T s = 0.4 and T w = 0.4 a,Q. 
These peak values correspond to the resonant photon energies. As seen from the figure, 
while the photoionization has very large value at Ri = 0.4 a^, it decreases too sharp with 
increasing Ri and reaches the its smallest value at Ri = 0.8 Oq. When the core radius 
is Ri = 0.9 <2q, it increases instantly and remains stable with further increasing of the 
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FIG. 7. The photoionization cross section of the hydrogenic donor impurity vs. core radius at 
constant T s and T w . 

R\. The physical reason of this treatment is as follows: at very small R\ values, the wave 
functions of both ground and excited states are localized in the well region as similar to 
Fig. H] and hence the overlapping of the wave functions becomes larger. As a result of 
this, the photoionization cross section is also larger. In this situation, the photoionization 
takes place in the well region and this may be important for a photodetector application 
which is sensitive to higher photon energies (blue shift) relatively, because the energy levels 
become larger, in smaller core radii. The increase in the core radius causes a decrease in 
the overlapping of the wave functions because the impurity wave function is confined in 
the core region while the excited states wave function is still confined in the well region. 
Therefore, the photoionization cross section has small values. With further increasing of the 
core radius, both ground and excited states wave functions are localized in the core region 
and hence the overlapping of the wave functions becomes large and as depending on this, 
the cross section is also to be maximum. In this case, the photoionization occurs in the core 
region. This situation may be important to produce a photodetector which is sensitive to 
lower energies (red shift). 

Figure M shows the variation of the photoionization cross section as a function of the shell 
thickness, T s , for two different values of the core radius, R\, and the same well thickness, 
T w = 0.4 Oq. When we look at the figure, we see that the cross section increases with in- 
creasing shell thickness for R± = 0.3 a^. On the other hand, for Ri = 0.5 <2q, in contrast with 
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FIG. 8. (Color online) The variation of the photoionization cross section of the impurity depending 
on shell thickness for two different core radii and constant T w . 

the previous case, the cross section exhibits completely different behavior and it decreases 
rapidly with increasing T s . Similar treatments are also observed for i? x > 0.5 values. 
When Ri = 0.3 a^, the wave functions of ground and excited states are confined in the well 
region and the finding probability in the core region vanishes almost for large values of the 
shell thickness. Therefore, the overlap between the wave functions becomes much larger in 
the well region as similar to the left panel of Fig. |H This is the reason of an increase in the 
photoionization cross section which happens in the well region. When Ri = 0.5 Oq, ground 
state wave function is confined in the core region because of the attractive coulomb potential 
of the impurity. At small values of T s , the tunneling of the excited state wave function to 
the core region is much more. In this case, the overlapping of the wave functions becomes 
larger and hence the photoionization cross section becomes greater. When the barrier thick- 
ness increases, the overlapping of the wave functions decreases because of the localizing of 
excited state wave function in the well region and therefore the cross section becomes less. 
We conclude that the sensitivity of a device can be also tuned with the shell thickness of 
MLQD to a narrower or larger region. 

The variation of the photoionization cross section with the well width, T w , is seen in Fig. 
[9] for R\ = 0.8 dp and two different shell thickness, T s = 0.2 and 0.4 a^. As seen from the 
figure, in small well widths, the photoionization cross section remains a bit stable and then 
it decreases rapidly with increasing well widths and reaches the minimum values at about 
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FIG. 9. (Color online) The variation of the photoionization cross section of the impurity as a 
function of well width for two different shell thickness and constant core radius. 

T w = 0.6 Oq. The minimum value of the photoionization for T s = 0.4 Oq becomes smaller 
than that for T s = 0.2 a^. The cross section increases again for both cases with further 
increasing of the T w . We can explain this behavior as follows: In small T w values, namely 
narrow well region case, the wave functions of ground and excited states are localized in 
the core region and therefore the overlapping of them becomes greater. The photoionization 
occurs in the core region. The excited state wave function starts to tunnel through the 
well region with increasing well width and it is completely localized in this region when the 
well width is about T w = 0.6 as seen from Fig. [TUJ In this case, the overlapping of the 
wave functions is minimum and hence, the photoionization cross section is very small. This 
decreasing is much more for T s = 0.4 Oq, because the tunneling of excited state wave function 
to the core region is smaller as shown at bottom panel of Fig. [TUJ By further increase of the 
well width, the ground state wave function also begins to tunnel to the well region and as 
a result of this, the overlapping of the wave functions increases again. When T w > 1.0 a^, 
both of the wave functions are localized in the well region and the cross section is to be 
higher. The photoionization cross section occurs in the well region anymore. 

The investigation of the photoionization cross section as depending on the photon energy 
is another important process. Figure [TT] shows the variation of the photoionization cross 
section of the neutral donor impurity as a function of the normalized photon energy hw/(Ef — 
Ei) for different core radii, R\, and shell thickness, T s . As seen from the figure, the maximum 
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FIG. 10. (Color online) Localization and energy cases for ground state of Z=l (black lines) and 
excited state of Z=0 (red lines). The left panel for T s = 0.2 clq and the right panel for T s = 0.4 clq. 



photoionization corresponds to the threshold photon energy for all R± and T s values and it 
decreases with increasing photon energy. These general trends are quite similar to those in 
single QDs.-^^ 1 ^ Also, the cross sections go towards lower values with increasing R\ and 
decreasing T s . When R\ is small, both ground and excited state wave functions are confined 
in the well region which results in a strong overlap between the wave functions and therefore 
in a larger photoionization cross section. The increase in the core radius, which corresponds 
to the slimming of the shell wall, leads to a possibility of tunneling of ground state wave 
function through the core region due to the attractive Coulomb potential of the impurity. 
So this possibility reduces the value of dipole matrix element and also as depending on this, 
the cross section becomes smaller with increasing normalized photon energy. 

In Fig. [121 variation of the photoionization cross section with normalized photon energy 
is shown for different well widths at R\ = 0.6 <2q and T s = 0.2<2q. As seen from the figure, the 
maximum photoionization cross section corresponds to the threshold energy and decreases 
with increasing photon energies for all well widths. When the well width increases, the 
photoionization cross section increases and this increase is more obvious for T w > 0.4 a^. 
This is originated from the overlapping of the wave functions in the well region. That is, 
the photoionization happens in the well region. 

As seen from Figs. (TT] and [12], the maximum absorption occurs at the photoionization 
threshold and the photoionization cross section decreases monotonically from the maximum 
value with increasing photon energy like the true hydrogenic model. This transition requires 
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FIG. 11. (Color online) Photoionization cross section of the hydrogenic donor impurity as a function 
of normalized photon energy for different R\, and T s and constant T w . 




FIG. 12. (Color online) Photoionization cross section of the hydrogenic donor impurity as a function 
of normalized photon energy for different T w and constant R\, and T s . 



higher photon energies than those in quantum wells and quantum wires.^ 1 ^ The reason for 
this behavior is that the electron is completely confined in all spatial dimensions in QD 
structures. This interesting result, which seems to be the characteristic and the signature 
of the QDs, is not observed in quantum well and quantum wire structures.-^>2& 
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CONCLUSION 



In this study, a detailed examination of the electronic structure and photoionization 
cross section of a donor impurity have been performed for a multi-layered spherical QD. In 
calculations, we have used shooting method with finite difference techniques and effective 
mass approximation for determining energy eigenvalues and their wave functions. First, we 
have calculated the donor binding energies as a function of the layer thickness. We have 
found that the donor binding energy strongly depends on the layer thickness. Then, we have 
computed photoionization cross section of the donor impurity as a function of layer thickness 
and normalized photon energies. We have showed that the photoionization cross section 
can be controlled with adjusting the core radius or layer widths. These properties can be 
useful for producing some device applications such as quantum dot infrared photodetectors. 
Although there is no experimental or theoretical reports in the literature to compare our 
results, general trend of our results are in a good agreement with previous studies of single 
quantum dots. We believe that this study will be rather useful and will contribute to the 
understanding of the photoionization cross section of MLQDs. We hope that this study will 
also stimulate both experimental and theoretical investigations of the photoionization cross 
section of MLQDs. 
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